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(54) Adaptive colour matching method and apparatus 



(57) Color matching method in which the color 
appearance of a color image is adapted to be coincident 
between color image devices having significantly differ- 
ent reference whites with a small processing volume. 
From correlated color temperature of an originating ref- 
erence white and a target reference white, spectral 
power distribution characteristics of an originating and a 
target color space are restored. From tristimulus values 
of the originating reference white, spectral power distri- 
bution characteristics of the originating color space (tar- 
get color space) and the human color matching 
functions, first (and second) white surface reflectances 
are restored. Through interpolation of the first and sec- 
ond white surface reflectances, an adaptation white sur- 
face reflectance is obtained. The ratio of the first white 
surface reflectance to the adaptation white surface 
reflectance is calculated to obtain a spectral chromatic 
adaptation ratio. The surface reflectance of the optional 
color in the originating color space is restored from the 
tristimulus values of an optional color, spectral power 
distribution characteristics of the originating color space 
and the human color matching functions. The surface 
reflectance of the optional color is multiplied by the 
spectral chromatic adaptation ratio to obtain an adap- 
tive surface reflectance of the optional color, and tristim- 
ulus values of the color in the target color space is 
obtained from the adaptive surface reflectance of the 
optional color, spectral power distribution characteris- 
tics of the target color space and the human color 
matching functions. 
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Description 

FIELD OF THE INVENTION 

[0001] This invention relates to a color reproduction technique for color images among different color image 
devices, such as a color scanner, a color monitor, a color printer etc. More particularly, it relates to a color matching 
method and apparatus required for converting an optional color in an originating color space to that of a target color 
space, in meeting with the chromatic adaptation state of the human visual system, as the color appearance correspond- 
ence is maintained between different color systems having different white references. 

BACKGROND OF THE INVENTION 

[0002] When a color image displayed on a monitor such as a cathode-ray tube (CRT) monitor for color reproduction 
between different color image devices is to be output to a color printer, color reproduction frequently presents a prob- 
lem. In general, the correlated color temperature of the reference white, set as a default reference in most CRT moni- 
tors, is set to approximately 9000K. In the case of a printer, on the other hand, the color when illuminated with D50 
(correlated color temperature: 5003K), as the standard light source for color evaluation for printing, is evaluated. In 
short, for color reproduction between monitor and printer, the color displayed on a monitor set to the standard white of 
9000K needs to be compared to the color of the printer output illuminated with the standard light source of 5000K. How- 
ever, even insofar as only the white is concerned, the appearance of the color with the correlated color temperature of 
9000 K perceived by us differs from that with the correlated color temperature of 5000 K. 

[0003] For sensibly realizing this difference between the two whites, two monitors are provided, with one of them 
being set to a standard white with a correlated color temperature of 9000K, with the remaining monitor being set to a 
standard white with a correlated color temperature of 5000K. If the respective whites are displayed on the two monitors 
and compared to each other, the difference in the color appearance can be perceived easily. If, in the state where the 
user's eye is completely adapted to the monitor of the standard white with the correlated color temperature of 9000K, 
the user sees the white displayed on the monitor of the standard white with the correlated color temperature of 5000K. 
the latter is perceived as being the yellowish color clearly different from the white. The reverse is also true, that is. if the 
user's eye is completely adapted to the monitor of the standard white with the correlated color temperature of 5000K, 
the user sees the white with the correlated color temperature of 9000K as being palish (bluish) white. 
[0004] However, if, in an environment that can control the ambient light, such as in a darkroom and, with a partition 
placed between both eyes, the monitor of 9000K and that of 5000K continue to be viewed with left and right eyes, 
respectively, the color appearance of the two whites gradually becomes similar (approach each other) due to progres- 
sive chromatic adaptation of both eyes. That is, since the color appearances of the two whites differ due to the state of 
our chromatic adaptation, such color conversion which takes our chromatic adaptation into account is indispensable in 
order to realize color coincidence between color image devices having difference reference whites. 
[0005] As a color conversion method, which takes the human chromatic adaptation into account, the von Kries 
model is well-known. This conversion method is shown in Fig. 10. According to this method, the chromatic adaptation 
is carried out based on the change in the RGB spectral sensitivity in the human visual system such that the RGB spec- 
tral sensitivity is upon change in illumination, changed in its sensitivity balance without changes in shape of the spectral 
curves so as to bring the two whites into coincidence. Assume RGB values of the illumination 1 be (R0, GO, B0), the 
RGB values of an object (article) in illumination 1 be (R, G, B), and RGB values of illumination 2 be (R0\ G0' f B0'), with 
RGB values of the same object under illumination 2 being (R\ G\ B* ), the tri-color sensual quantities of the visual sys- 
tem of the object color are expressed by: R/R0, G/G0, B/B0, R' /R0' , G* /GO' , B '/BO* . In order for the color appearance 
of the object to be coincident under the illumination 1 and under the illumination 2, it suffices if the above-mentioned tri- 
color sensual quantities are coincident, as shown below: 
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[0006] The RGB values can be obtained by linear transformation of tristimulus values X, Y, 2: 
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[0007] It is noted that RO, GO, BO. R0* , G0\ BO* in the above formula (1) may be obtained by substituting tristimulus 
values (XO. YO. ZO). (X0\ Y0\ ZO' ) of the illumination 1 and the illumination 2 shown in Fig. 1 1 into formula (3). 
[0008] In the XYZ->RGB transformation matrix M of formula (3). Pitt's matrix or Estevez' matrix may be used. 
[0009] By substituting equation (3) into equation (2), the following von Kries chromatic adaptation prediction for- 
w mula is obtained: 
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[001 0] In the above formula (4) , 
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[0011] It is seen from above that the color on the reproducing side corresponding to tristimulus values (X, Y, Z) of 
25 the input color on the original side in an observing booth shown in Fig. 1 1 can be calculated from the von Kries chro- 
matic adaptation prediction formula (4). 

[0012] Nayatani et al proposed a chromatic adaptation model combining the linear process in which the von Kries 
chromatic adaptation model is valid and the non-linear process in which an exponent varies responsive to the adapta- 
tion level. This chromatic adaptation model was recommended by CIE in 1986 as a chromatic adaptation prediction for- 
30 mula. This chromatic adaptation prediction formula is hereinafter explained. 

[0013] First, tristimulus values X, Y t Z of the input color on the original side in Fig. 11 are transformed, by using 
Estevez' matrix, into basic stimulus values RGB in accordance with the following formula (6): 
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[0014] From the basic stimulus values RGB of the input color, basic stimulus values R*. G\ B* of the corresponding 
40 colors are calculated. 

/?'=(100p 0 ^+1){(H + 1)/(100p 0 ^ + 1)} Pr -1 (7) 
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[001 5] For the corresponding colors R\ G\ B\ the formula (6) is inverse transformed (as shown by the following for- 
mula (8)) to calculate tristimulus values X\ Y\ Z' of the corresponding color: 
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[001 6] In the above formula (7), po, po' denote reflectances of the surrounding environments on the original and 
reproducing sides, respectively, with 0.2^po^0.1 and 0.2=ipo^1 .0, and n, and £ are as indicated by the following for- 
mulas: 

5 4=(0.48105x+0.78841y-0.08081)/y <9) 

ri =(-0.27200 x+1 .1 1962y+0.04570)/y 
£=0.91 822(1 -x-y)/y 

10 

where x, y denote chromaticity of the illumination 1 . £\ V , and £ are calculated in a similar manner using the chroma- 
ticity of the illumination 2 in x and y. 

Pr=f(R 0 yf{R Q <) (10) 

15 

P9=HG 0 )/f(G Q ') 



Pb=g(B Q )/g(B 0 ') 

20 where functions fO and g0 are as follows: 

/(*H6.469+6.362x° 4495 )/(6.469+x° 4495 ) (1 1) 

S(x)=(8.41 4+8.091 x°' 5128 )/(8.41 4+x 0 - 5128 ) (12) 

25 

[001 7] Nayatani et al improved the above-mentioned CIE1 986 chromatic adaptation prediction formula in 1 990 and 
in 1994 and, in 1995, completed a color appearance model (Nayatani model) in a publication (Y. Nayatani et al: Light- 
ness Dependency of Chroma Scales of a Nonlinear Color-Appearance Model and its Latest Formulation. Color Res. 
Appl., Vol. 20, No. 3, pp. 156-167, June 1995). 
30 Further the correction method for taking into account of incomplete chromatic adaptation state of the human visual sys- 
tem was added in 1996 and 1997. 

[0018] Fairchild has proposed RLAB as a color appearance model. The calculating method of the present model is 
explained. First, basic stimulus values LMS are found fromtristimulus values XY2 in the original side in accordance with 
the following formula (13): 

35 
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[0019] A transformation matrix A for chromatic adaptation then is found. The calculating method is as follows: 
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P s + D(1.0-P s ) 
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[0020] The tristimulus values XYZ of the input color are transformed into corresponding stimulus values Xref, Yref 
and Zref under the reference observation conditions in accordance with the following formula (24): 



Yref 
l Z ret. 



X 

= RAM Y 
Z. 



(24) 



[0021] In the above formulas, Yn is an absolute adaptation luminance of the illumination (cd/m 2 ), and Ln, Mn and 
Sn are the basic stimulus values of the illumination. The constant D in formulas (15) to (17) is set to 1 .0, 0.0 and to an 
intermediate value for a hard copy image, monitor display and for a slide, respectively. The constant matrix R in formula 
ao (24) is as follows: 



45 



so 
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-1.9569 -1.1882 0.2313H 
0.3612 0.6388 0.0 
L 0.0 0.0 1.0 



[0022] The RLAB values then are calculated using the reference tristimulus values Xref, Yref and Zref. 

L R =W(Y nf ) a 

a R =430[(X ret ) c -(Y rel )°] 

b R ^70[(Y ref )^{Z ref n 



(25) 



(26) 
(27) 
(28) 



where a assumes the following values, responsive to the observing environment: 
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a=1/2.3 for an average ambient environment; 
o=1/2.9 for a dim ambient environment; and 
a=1/3.5 for a dark ambient environment. 

5 [0023] As tor the corresponding colors X'Y'Z' under the illumination 2, reverse calculations of the formulas (26) to 
(28) are carried out from the RLAB values corresponding to the input color of the illumination 1 to compute the reference 
tristimulus values Xref , Yref and Zref for the illumination 2. 

[0024] The corresponding colors X* Y' T in the illumination 2 can be calculated by substitution into the right side of 
the equation (29) which is the reverse calculation of the formula (24). It is noted however that the matrix calculated for 
70 the illumination 2 can be used as the matrix A of the formula (29). 

[0025] As other color appearance models, a Hunt model by Hunt, LLAB by Luo and CIECAM97s recommended by 
CIE in autumn of 1997, are attracting attention as models of significance, and evaluation is underway on the color 
appearance models among cross media. 

is SUMMARY OF THE DISCLOSURE 

[0026] As for the color appearance coincidence of the color images between color image devices with different ref- 
erence whites, such as between a monitor and a printer, there is a problem of color appearance to which the human 
chromatic adaptation as a human visual function is intricately pertinent. 

20 [0027] Many of the currently proposed color appearance models perform modeling of various color appearance 
phenomena of the human color perception through a complex calculation formula. In particular, Hunt model and Naya- 
tani model require very complicated calculations in order to represent the various color appearance phenomena. It is, 
however, open to doubt whether all the appearance phenomena considered in these models are really necessary in the 
practical utilization of color management on the computer system. That is, if the totality of appearance phenomena 

25 would be considered, the proposed chromatic adaptation models cannot be said to be practical in view of correspond- 
ingly increased computational costs. 

[0028] On the other hand, the von Kries chromatic adaptation model, which changes the ratio of the spectral sen- 
sitivity on the human retina level, responsive to the illumination color, without changing its shape of spectral curves, is 
suited to the color management function on the computer system, because of the smaller number of parameters 
30 required for calculation and the low computing costs. However, as for the problem of color appearance between a mon- 
itor and a printer having acutely different reference whites, it is not possible with the von Kries model to realize an opti- 
mum color appearance. 

[0029] There is a demand for a model which gives practically coincident color appearance of the color images and 
entails low computing costs in realizing the color matching of color images between different color image devices on a 
35 computer system. 

[0030] In view of the above-described status of the art, it is an object of the present invention to provide a method 
and apparatus for color matching according to which coincidence of color appearance of color images between color 
image devices having significantly different reference whites can be achieved with a small computational volume taking 
into account the chromatic adaptation state of the human visual system. 
40 [0031 ] According to a first aspect of the present invention, there is provided a color matching method for converting 
an optional color in an originating color space to a color in a target color space, keeping the color appearance. The 
method comprises the steps of: 

(a) restoring spectral power distribution characteristics of an originating color space from a correlated color temper- 
45 ature of an originating reference white which is the reference white in said originating color space, and restoring 

spectral power distribution characteristics of a target (objective) color space from a correlated color temperature of 
a target (objective) reference white which is the reference white in the target color space; 

(b) restoring, as a surface reflectance of the originating reference white, a first white surface reflectance from tris- 
timulus values of the originating reference white, spectral power distribution characteristics of the originating color 

so space and color matching functions of the human visual system (termed as " human color matching functions" 
hereinafter), and also restoring a second white surface reflectance from tristimulus values of the originating refer- 
ence white, spectral power distribution characteristics (i.e.. spectra) of the target color space and the human color 
matching functions; 

(c) obtaining an adaptive white surface reflectance thorough interpolation of the first and second white surface 
55 reflectances; 

(d) calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance to obtain a spec- 
tral chromatic adaptation ratio; 

(e) restoring a surface reflectance of the optional color in the originating color space by employing tristimulus values 
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of the optional color, spectral power distribution characteristics of the originating color space and the human color 
matching functions: 

(0 multiplying the surface reflectance of the optional color with the spectral chromatic adaptation ratio to obtain an 
adaptive surface reflectance of the optional color; and 
5 (g) finding tristimulus values of the color in the target color space from the adaptive surface reflectance of the 

optional color, the spectral power distribution characteristics of the target color space and the human color match- 
ing functions. 

[0032] With the color matching method of the present invention, in the step (c) of obtaining adaptation white surface 
to reflectance from the interpolation of the first and second white surface reflectances, an incomplete chromatic adapta- 
tion coefficient is changed responsive to the state of incomplete chromatic adaptation in the human visual system, the 
degree of interpolation of the first and second white surface reflectances is freely changed to obtain the adaptation 
white surface reflectance. 

[0033] The color matching method of the present invention further comprises the step of calculating an optimum 
is incomplete chromatic adaptation coefficient from combination of the reference white of the originating color space and 
the reference white of the target color space, the incomplete chromatic adaptation coefficient obtained from the opti- 
mum incomplete chromatic adaptation coefficient calculating step being used. 

[0034] According to a second aspect of the present invention, a color matching apparatus is configured for convert- 
ing an optional color in an originating color image device to a color in a target color image device, keeping the color 
so appearance. The color matching apparatus comprises: 

(a) means for restoring spectral power distribution characteristics of an originating color space from a correlated 
color temperature of an originating reference white which is the reference white of the originating color image 
device, and restoring spectral power distribution characteristics of the target color space from a correlated color 

25 temperature of a target reference white which is the reference white of the target color image device, 

(b) means for restoring, as a surface reflectance of the originating reference white, a first white surface reflectance 
from tristimulus values of the originating reference white, the spectral power distribution characteristics of the orig- 
inating color space and human cdor matching functions, and also restoring a second white surface reflectance 
from tristimulus values of the originating reference white, spectral power distribution characteristics of the target 

30 color space and the human color matching functions. 

(c) means for obtaining an adaptive white surface reflectance through interpolation of the first and second white 
surface reflectances, 

(d) means for calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance to 
obtain a spectral chromatic adaptation ratio. 

35 (e) means for restoring a surface reflectance of the optional color in the originating color image device by employing 
tristimulus values of the optional color, spectral power distribution characteristics of the originating color space and 
the human color matching functions. 

(f) means for multiplying the surface reflectance of the optional color with the spectral chromatic adaptation ratio to 
obtain an adaptive surface reflectance of the optional color, 
40 (g) means for finding tristimulus values of the color in the target color image device from the adaptive surface 
reflectance of the optional color, the spectral power distribution characteristics of the target color space and the 
human color matching functions. 

[0035] In the color matching apparatus of the present invention, an incomplete chromatic adaptation coefficient is 
45 changed in the means for obtaining adaptive white surface reflectance from the interpolation of the first and second 
white surface reflectances, responsive to the state of incomplete chromatic adaptation in the human visual system, so 
that the degree of interpolation of the first and second white surface reflectances freely changed to obtain the adaptive 
white surface reflectance. 

[0036] The color matching apparatus of the present invention also includes means for calculating an optimum 
so incomplete chromatic adaptation coefficient from combination of the reference white of the originating color image 
device and the reference white of the target color image device, the incomplete chromatic adaptation coefficient 
obtained from the optimum incomplete chromatic adaptation coefficient calculating means being then used. 
[0037] The color matching apparatus also includes means for computing the correlated color temperature of the 
reference white from the chromaticity of the reference white of the color image device. 

55 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0038] 

5 Fig. 1 is a block diagram showing the structure of a first embodiment of the color matching device according to the 

present invention. 

Fig. 2 is a block diagram showing the structure of a spectral chromatic adaptation ratio computing means in the first 
embodiment of the present invention. 

Fig. 3 is a block diagram showing the structure of a second embodiment of the color matching device according to 
ro the present invention. 

Fig. 4 is a block diagram showing the structure of a third embodiment of the color matching device according to the 
present invention. 

Fig. 5 is a block diagram showing an automatic calculating device for incomplete chromatic adaptation coefficient in 
the third embodiment of the present invention. 
75 Fig.6 is a block diagram showing the structure of a fourth embodiment of the color matching device according to 

the present invention. 

Fig.7 is a block diagram showing the structure of a fifth embodiment of the color matching device according to the 
present invention. 

Fig.8 is a block diagram showing an illustrative structure of the hardware configuration for realizing a color matching 
20 device according to the present invention. 

Fig.9 is a flow diagram showing processing flow of the color matching device according to the present invention. 
Fig. 10 illustrates a von Kries chromatic adaptation model. 

Fig.11 illustrates the state of color comparison under different illuminating environments. 

25 PREFERRED EMBODIMENTS OF THE INVENTION 

[0039] Preferred embodiments of the present invention are explained. In its preferred form, the color matching 
apparatus of the present invention includes hypothetical spectral power distribution characteristics calculating means 
(100 of Fig.1) for restoring spectral power distribution characteristics (3 of Fig.1) of the originating color space from the 
30 correlated color temperature (1 of Fig.1) of the originating reference white as a reference white in the originating color 
image device (201 of Fig.1), and for restoring spectral power distribution characteristics (4 of Fig.1) in the target color 
space from the correlated color temperature (2 of Fig. 1 ) of the target reference white as the reference white in the target 
color image device (202 of Fig.1). The color matching apparatus of the present invention also includes hypothetical sur- 
face reflectance calculating means (10V of Fig.2) for restoring, as a surface reflectance of the originating reference 
35 white, a first white surface reflectance (1 1 of Fig.2) from tristimulus values of the originating reference white (10 of 
Fig.1), spectral power distribution characteristics of the originating color space (3 of Fig.1) and human color matching 
functions, and for restoring a second white surface reflectance (12 of Fig.2) from tristimulus values of the originating ref- 
erence white (10 of Fig.1), spectral power distribution characteristics of the target color space (4 of Fig.1) and human 
color matching functions. The color matching apparatus of the present invention also includes incomplete absolute 
40 white adaptation calculating means (107 of Fig.2) for obtaining adaptive white surface reflectance (14 of Fig.2) based 
on interpolation of the first and second white surface reflectances, and ratio computing means (108 of Fig.2) for obtain- 
ing the spectral chromatic adaptation ratio (7 of Figs. 1 and 2) based on calculations of the ratio of the first white surface 
reflectance and adaptive white surface reflectance. The color matching apparatus of the present invention also includes 
hypothetical reflectance calculating means (101 of Fig.1) for restoring surface reflectance of an optional color in an orig- 
45 inating color image device (6 in Fig.1) by using tristimulus values of the optional color (5 of Fig.1), spectral power distri- 
bution characteristics (3 of Fig.1) of the originating color space and human color matching functions. The color 
matching apparatus of the present invention also includes incomplete chromatic adaptation processing means (104 of 
Fig.1) for obtaining the adaptive surface reflectance (8 of Fig.1) of the optional color by multiplying the surface reflect- 
ance of the optional color (6 in Fig.1) with a spectral chromatic adaptation ratio (7 of Fig.1) and tristimulus values corn- 
so puting means (106 of Fig.1) for finding the tristimulus values of the color in the target color image device from the 
adaptive surface reflectance of the optional color, spectral power distribution characteristics of the target color space 
and human color matching functions. The color matching apparatus of the present invention also includes XYZ-RGB 
converting means (106 of Fig.1) for converting the tristimulus values of the target color image device into device color 
data of the target color image device. 
55 [0040] In its alternative form, the color matching apparatus of the present invention includes means for calculating 
an optimum incomplete chromatic adaptation coefficient (206 of Fig.5) for finding an optimum incomplete chromatic 
adaptation coefficient (13 of Fig.5) from combination of correlated color temperature of the reference white (1 of Fig.5) 
of the originating color image device and correlated color temperature (2 of Fig.5) of the target reference white of the 
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target color image device, the incomplete chromatic adaptation coefficient obtained from the optimum incomplete chro- 
matic adaptation coefficient calculating means being used. 

[0041] The color matching apparatus of the present invention also includes calculating means (109 of Fig.6) for cal- 
culating the correlated color temperature of the reference white from the chromaticfty of the reference white of the color 
5 image device. 

[0042] The color matching method of the present invention includes a step of restoring spectral power distribution 
characteristics of an originating color space from a correlated cdor temperature of an originating reference white which 
is the reference white of the originating color space and restoring spectral power cfistribution characteristics of a target 
color space from a correlated color temperature of a target reference white which is the reference white of the target 
w color space (S1 and S2 of Fig.9); 

a step of restoring, as a surface reflectance of the originating reference white, a first white surface reflectance from 
tristimulus values of the originating reference white, the spectral power distribution characteristics of the originating 
color space and human color matching functions, and also restoring a second white surface reflectance from tris- 
15 timulus values of the originating reference white, spectral power distribution characteristics of the target color space 
and human color matching functions (S3 and S4 in Fig. 9); 

a step of obtaining an adaptive white surface reflectance through interpolation of the first and second white surface 
reflectances (S5 of Fig.9); 

a step of calculating the ratio of the first white surface reflectance to the adaptive white surface reflectance to obtain 
20 a spectral chromatic adaptation ratio (S6 of Fig.9); 

a step of restoring a surface reflectance of the optional color in the originating color space by employing tristimulus 
values of the optional color, spectral power distribution characteristics of the originating color space and human 
color matching functions (spectral chromatic adaptation ratio) (S7 of Fig.9); 

a step of multiplying the surface reflectance of the optional color with the spectral chromatic adaptation ratio to 
25 obtain an adaptive surface reflectance of the optional color (S8 of Fig.9); and 

a step of finding tristimulus values of the color in the target color space from the adaptive surface reflectance of the 
optional color, the spectral power distribution characteristics of the target color space and the human color match- 
ing function (S9 of Fig.9). 

30 [0043] In the present embodiment of the present invention, the processing of steps S1 to S9 may be carried out by 
program control executed on a computer. The present invention can be carried out by reading out a program from a 
recording medium having recorded the program thereon and executing the read-out program. This program is deemed 
a computer readable program product which may be also carried/contained/transmitted by any medium available which 
includes static and dynamic media, e.g., recording medium, physical medium such as carrier wave of electromagnetic 

35 and/or optical nature. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

[0044] Referring to the drawings, present embodiments of the present invention will be explained in detail. First, an 
40 embodiment of the color matching method according to the present invention is explained. 

[0045] The color matching in case of two color image devices, that is a color image device on the original side and 
another one on the reproducing side, have different reference whites, is explained. It is assumed that the correlated 
color temperature of the reference whites of the original side and the reproducing side color image device and tristim- 
ulus values XYZ are given from the outset. 
45 [0046] First, hypothetical spectral power distribution characteristics lv1(X), equivalent to the color of illumination 1 
in a color image scene reproduced by a color image device on the original side, and hypothetical spectral power distri- 
bution characteristics lv2(X), equivalent to the color of illumination 2 in a color image scene reproduced by a color image 
device on the reproducing side, are calculated based on a correlated color temperature of a reference white light of 
each color image device. As the method of calculating the hypothetical spectral power distribution characteristics, 
so equivalent to the illumination color in the color image scene, the technique disclosed in the Japanese Patent Application 
H-9-047334 (Kokai Publication J P -A- 10-229499) by the present inventor, may be used, the entire disclosure of which is 
incorporated herein by reference thereto. 

[0047] Given an absolute white as an input color on the original side, its tristimulus values (Xw, Yw. Zw) coincide 
with tristimulus values (X0. Y0. ZD) of the illumination 1 in the image scene. 

55 

(Xw. Yw, Zw) = (X0. Y0, Z0) 

[0048] From the tristimulus values (X0. Y0, Z0) of the absolute white and the hypothetical spectral power distribu- 
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tion characteristics Iv1 (X) calculated from the correlated color temperature of the illumination 1 , It is possible to calculate 
the hypothetical surface reflectance Rwv(A.) of the absolute white on the original side. For restoring the hypothetical sur- 
face reflectance of an object from the tristimulus values of the object and hypothetical spectral power distribution char- 
acteristics of the illumination, a restoration method described in the Japanese patent Application H-9-047334 may be 
5 used. 

[0049] It is now envisaged to reproduce the tristimulus values (Xw, Yw, Zw) of the absolute white on the original 
side on a reference side color image device having the illumination 2 in the scene with the identical tristimulus values 
(Xw, Yw. Zw). From the tristimulus values (Xw, Yw, Zw) of the absolute white on the original side and the hypothetical 
spectral power distribution characteristics lv2(X.) calculated from the correlated color temperature of the illumination 2, 
10 the hypothetical surface reflectance R'wv(X) of the absolute white on the reference side is calculated using a method 
similarly described in the Japanese Patent Application H-9-047334 (JP-A-1 0-229499). 

[0050] It is seen from above that, as to the absolute white, two hypothetical surface reflectances, namely Rwv(A.) on 
the original side color image device and R'wvfA.) on the reference side color image device have now been obtained. 
[0051] For coping with the situation in which the color perception of the viewer has been incompletely adapted 
15 between the original and reference side color image devices, the incomplete chromatic adaptation coefficient adj (real 
number such that 0 ^ adj = 1) is defined. 

[0052] The hypothetical surface reflectance Rwv ad (A,). for which the color appearance of the absolute white of the 
two color image devices coincide in the state of incomplete chromatic adaptation of a viewer, is calculated, using the 
calculated hypothetical surface reflectances Rwv(X) and Rwv'(X.) of the absolute white, in accordance with the following 
20 formula (29): 

R\w ad {X)=adjxRwv(X)+(-\-adj)xRwv'(X) (29) 

[0053] As for colors other than the absolute white, the following explanation is given. From the hypothetical surface 
25 reflectance Rwv ad (A.) of the absolute white in the state of incomplete chromatic adaptation obtained by the formula (29) 
and the hypothetical surface reflectance Rwv(A.) of the absolute white under the illumination 1 , a specific coefficient is 
calculated for each wavelength in the visible light range and represented as ref^f (X) as the spectral chromatic adap- 
tation ratio. The spectral chromatic adaptation ratio ref COGf (X) is represented by the following formula (30): 

30 ref wef (X)^Rwv ad (X)IRwv(X) (30) 

[0054] From the hypothetical spectral power distribution characteristics Iv1 (X), calculated from the correlated color 
temperature of the illumination 1 , and the tristimulus values XYZ of an optional color in the color image on the original 
side, the hypothetical surface reflectance Rv(X) equivalent for the object color is calculated, using the technique stated 
35 in, for example, Japanese Patent Application H-9-047334 (JP-A-1 0-229499). 

[0055] The hypothetical surface reflectance Rv(X) is multiplied with the spectral chromatic adaptation ratio ref coef 
(X) of the formula 30, from one wavelength to another. That is. the hypothetical surface reflectance Rv ad (X) of the 
optional input color in the incomplete chromatic adaptation state is given by the following formula (31): 

40 Rv ad (X)= Rv(X) x ref coe/ (X) (31 ) 

[0056] Tristimulus values X\ Y\ Z' of a corresponding color on the reproducing side color image device, which color 
corresponds to an input color in a color image, reproduced on the original side color image device, can be calculated 
by the following formula (32): 

45 

XW/*(X)ffl' acr (X)j(X)cA. (32) 

Y'-UygWRv^XMXXtk 

50 ^Hv2MRv ad (X)z(X)dX 

exploiting the hypothetical surface reflectance Rv ad (A.) of the object in the state of incomplete chromatic adaptation and 
the calculated hypothetical spectral power distribution characteristics lv2(A.) of the illumination on the reproducing side. 
[0057] In the formula (32), x (X), y (X), and z {X) are known color matching functions. 
55 [0058] An embodiment of a color matching device according to the present invention is now explained. Fig.1 shows 
a block diagram showing the structure of a first embodiment of the color matching device according to the present 
invention. 

[0059] Referring to Fig.1, to a color matching device 200 are connected a color image display device 201 on the 
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original side for color reproduction with RGB values and a color image cfisplay device 202 on the reference side. The 
color characteristics of the respective color image display devices are assumed to be known from the outset The color 
matching device 200 converts an optional color reproduced on the original side color image display device 201 into a 
color on the reproducing side color image display device 202 in such a manner as to achieve coincident color appear- 
5 ance of the two colors. 

[0060] The color matching device 200 includes a hypothetical spectral power distrfoution characteristics calculating 
means 100, a hypothetical reflectance calculating means 101, an RGB-XYZ conversion means 102, a spectral chro- 
matic adaptation ratio calculating means 103. an incomplete chromatic adaptation processing means 104, tristimulus 
values calculating means 105 and an XYZ-RGB conversion means 106. An incomplete chromatic adaptation coefficient 

w 13 is entered as an input. 

[0061 ] The operation of the f irst embodiment of the color matching device of the present invention is now explained. 
[0062] The correlated color temperature 1 of the original reference white of the original side color image display 
device 201 is given as color characteristics of the color image display device. This conelated color temperature 1 of the 
original reference white is converted by the hypothetical spectral power distribution characteristics (spectra) calculating 

15 means 100 into hypothetical spectral power distribution characteristics (spectra) 3 of the originating color space corre- 
sponding to the illumination color of the original side. For this hypothetical spectral power distribution characteristics cal- 
culating means 100, the technique disclosed in for example the Japanese patent application H-9-047334 may, for 
example, be used. 

[0063] Similarly, a correlated color temperature 2 of the target reference white of the reproducing side color image 
20 display device 202 is also given as color characteristics, and is converted by the hypothetical spectral power distribution 
characteristics calculating means 100 into hypothetical spectral power distribution characteristics (spectra) 4 of the tar- 
get color space corresponding to the illumination color of the reproducing side. 

[0064] The RGB-XYZ conversion means 102 converts optional RGB values in the color space of the color image 
display device 201 of the original side into tristimulus values XYZ. Since the color characteristics of the original side 
25 color image display device 201 are known from the outset, the RGB values may be converted by any suitable known 
technique into the tristimulus values XYZ. In this manner. XYZ values 5 of the input color are calculated for an optional 
color of the color space reproduced by the color image display device 201 of the original side. 

[0065] From the hypothetical spectral power distribution characteristics 3 of the originating color space correspond- 
ing to the of the illumination color on the original side, calculated by the hypothetical spectral power distribution charac- 

30 teristics calculating means 1 00, and the input color XYZ values 5 from the color image reproduced on the original side 
color image display device 201. the hypothetical reflectance calculating means 101 calculates a hypothetical surface 
reflectance 6 of the input color corresponding to the object color not dependent on the illumination. For the hypothetical 
reflectance calculating means 101, the technique disclosed in, for example, the Japanese patent application H-9- 
047334 (JP-A-1 0-229499) etc. may be used. 

35 [0066] The spectral chromatic adaptation ratio calculating means 103 calculates a spectral chromatic adaptation 
ratio ref coef (X), used for realizing color appearance matching in the state of incomplete chromatic adaptation of the 
human color perception on the color matching method of the present invention. 

[0067] Fig.2 shows a block diagram showing the detailed structure of the spectral chromatic adaptation ratio calcu- 
lating means 103. Referring to Fig.2, the operation of the spectral chromatic adaptation ratio calculating means 103 is 
40 explained. 

[0068] Referring to Fig.2. the spectral chromatic adaptation ratio calculating means 103 calculates a first hypothet- 
ical white surface reflectance 11, in a hypothetical surface reflectance calculating means 10V, from tristimulus values 
XYZ 10 of the absolute white, reproduced on the color image display device 201 on the original side and hypothetical 
spectral power distribution characteristics 3 of the originating color space corresponding to the illumination color on the 

45 original side. This hypothetical surface reflectance calculating means 10V is the same as the hypothetical reflectance 
calculating means 101 explained in connection with the color matching device 200. Also, a second hypothetical white 
surface reflectance 12 is calculated from the tristimulus values XYZ 10 of the absolute white, reproduced on the color 
image display device 201 on the original side, and hypothetical spectral power distribution characteristics 4 of the target 
color space, corresponding to the color of the illumination on the reproducing side. 

so [0069] An incomplete chromatic adaptation calculating means for absolute white 107 calculates, by the formula 
(29). the hypothetical surface reflectance 14 of the absolute white in the incomplete chromatic adaptation state for 
which color appearance of the absolute white in the incomplete chromatic adaptation state of the viewer coincides 
between the original and reference side color image devices, with the first hypothetical white surface reflectance 1 1 , 
second hypothetical white surface reflectance 12 and an incomplete chromatic adaptation coefficient 13 as inputs. The 

55 incomplete chromatic adaptation coefficient 13 is used for coping with (or compensating for) the state of incomplete 
chromatic adaptation of the visual perception of the viewer between the original side color image display device and the 
reference side color image display device. 

[0070] A ratio computing means 108 calculates spectral chromatic adaptation ratio 7 from a hypothetical surface 
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reflectance 1 4 of the absolute white in the incomplete chromatic adaptation state and the first hypothetical white surface 
reflectance 1 1 . The computing method is as shown in the formula (30). 

[0071] Reference is again had to Fig.1 for reverting to the explanation of the color matching device 200. The incom- 
plete chromatic adaptation processing means 104 multiplies the hypothetical surface reflectance 6 of the input color on 
5 the original side by the above-mentioned spectral chromatic adaptation ratio 7, to output the hypothetical surface 
reflectance of the optional input color in the incomplete chromatic adaptation state, that is a hypothetical surface reflect- 
ance 8 of the corresponding color. 

[0072] The tristimulus value calculating means 1 05 computes X\ Y\ 2* of the corresponding color 9, in accordance 
with the formula (32), from the hypothetical surface reflectance 8 of the corresponding color, hypothetical spectral power 
10 distribution characteristics 4 corresponding to the illumination color on the reproducing side and from the color matching 
functions. 

[0073] Finally, the X\ Y\ Z of the corresponding color 9 are converted by the XYZ-RGB conversion means 1 06 into 
RGB value of the color image display device 202 on the reproducing side to output the resulting RGB values to the color 
image display device 202 on the reproducing side. 
15 [0074] Fig.3 shows the block diagram showing the structure of a second embodiment of the color matching appa- 
ratus according to the present invention. In Fig.3, the parts or components which are the same as those shown in Fig.1 
are denoted by the same reference numerals. 

[0075] Referring to Fig.3 f a color image input device 203 on the original side and a color image display device 202 
on the reference side are connected to a color matching device 204. The color image input device 203 on the original 
20 side is used for reading the color document to formulate color image data constituted by RGB values. Meanwhile, the 
color characteristics of the color image input device 203 on the original side and the color image display device 202 are 
assumed to be known from the outset. 

[0076] Referring to Fig.3, this color matching device 204 includes a correlated color temperature conputing means 
109, in addition to the structure provided in the color matching device 200 shown in Fig.1 . The correlated color temper- 
as ature computing means 109 finds the correlated color temperature 1 of the originating reference white based on the 
chromaticity of the reference white of the color image input device 203 on the input side. As the method for computing 
the corresponding correlated color temperature from the chromaticity of the reference white used in the correlated color 
temperature computing means 109, the technique described in, for example, the Japanese patent application H-9- 
047334 (J P-A- 10-229499) may. for example, be used. 
so [0077] The correlated color temperature 1 of the originating reference white, calculated by the correlated color tem- 
perature computing means 1 09, is input to the hypothetical spectral power distribution characteristics calculating means 
1 00. The processing similar to that performed in the color matching device 200 shown in Fig. 1 then is carried out so that 
the optional input color captured by the color image input device 203 on the original side is converted into the color of 
the color image display device 202 of the reference side and output to the color image display device 202 on the refer- 
35 enceside. 

[0078] Fig.4 shows a block diagram showing the structure of a third embodiment of the color matching apparatus 
of the present invention. To a color matching apparatus 205, there are connected the color image display device 201 of 
the original side and the color image display device 202 of the reference side. The color characteristics of the respective 
color image display devices are assumed to be known from the outset. Referring to Fig.4, this color matching apparatus 
40 205 is made up of the color matching device 200 of the previous first embodiment and an automatic calculating device 
for incomplete chromatic adaptation coefficient 206. 

[0079] The correlated color temperature 1 of the originating reference white and the correlated color temperature 2 
of the target (object) reference white are obtained from the color characteristics of the color image display device. The 
correlated color temperature 1 of the originating reference white and the correlated color temperature 2 of the target 

45 reference white are input to the automatic calculating device for incomplete chromatic adaptation coefficient 206. 
[0080] The automatic calculating device for incomplete chromatic adaptation coefficient 206 is explained. 
[0081] Fig.5 shows a block diagram showing an illustrative structure of the automatic calculating device for incom- 
plete chromatic adaptation coefficient 206. Referring to Fig.5, the automatic calculating device for incomplete chromatic 
adaptation coefficient 206 includes an optimum incomplete chromatic adaptation coefficient calculating means 1 10 and 

so an incomplete chromatic adaptation coefficient memory 15. The automatic calculating device for incomplete chromatic 
adaptation coefficient 206 calculates an incomplete chromatic adaptation coefficient 13 for coping with (or compensat- 
ing for) a state in which the color perception of the viewer has adapted itself incompletely between the reference white 
of the original side color image device and that of the reference side color image device. 

[0082] First, an optimum incomplete chromatic adaptation coefficient adj is obtained for each of the various combi- 
55 nations of the correlated color temperatures, such as 9000K-6500K or 9000K-5000K. and is stored in the incomplete 
chromatic adaptation coefficient memory 15. The optimum incomplete chromatic adaptation coefficient calculating 
means 110 captures, from the combinations of the correlated color temperature 1 of the input original reference white 
and the correlated color temperature 2 of the target reference white, an incomplete chromatic adaptation coefficient cor- 
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responding to the input combination from the incomplete chromatic adaptation coefficient memory 15, and outputs an 
incomplete chromatic adaptation coefficient 13. This incomplete chromatic adaptation coefficient 13 is entered as an 
input to the color matching device 200. 

[0083] The process of converting an optional color of the color image display device 201 of the original side to a 
5 color of the color image display device 202 of the reference side by the color matching device 200 is as explained in the 
first embodiment described previously. 

[0084] The color image input device 203 may also be connected to the color image device on the original side in 
the color matching apparatus 205. 

[0055] Fig.6 shows a block diagram showing a structure of a fourth embodiment of the color matching apparatus of 
io the present invention. Referring to Fig.6 t the color matching apparatus 207 is realized by connecting a color image input 
device 203 as an original side color image device in the color matching apparatus 205, adding a correlated color tem- 
perature calculating means 109 for calculating the correlated color temperature 1 of the originating reference white, and 
by exchanging a color matching device 204 for the color matching device 200. 

[0086] In the above-described embodiment, the color matching devices or apparatus 200, 205 and 206 have been 
is explained for the cases in which a color image input device is used on the original side and a color image display device 
is used on the reference side. Alternatively, a color image recording device for reproducing the color by multi-colored 
inks, such as cyan, magenta, yellow and black inks, can be connected to the reference side color image device. 
[0087] Fig.7 shows a block diagram showing a structure of a fifth embodiment of the color matching device of the 
present invention. In Fig.7, parts or components which are the same as or similar to those shown in Fig.1 are denoted 
20 by the same reference numerals. 

[0088] Referring to Fig.7, the original side color image display device 201 and the reference side color image print 
device 207 are connected to this color matching device 208. The reference side color image print device 207 is adapted 
to output a color image by multi-colored inks, such as cyan, magenta, yellow and black inks. It is assumed that color 
characteristics of the original side color image display device 201 and the reference side color image print device 207 
25 are known from the outset, and colorimetric data 16 of the reference side color image print device, obtained on meas- 
uring the relation between the amount of the multi-colored ink and the tristimulus values, are given. The color matching 
device 208 is modified from the color matching device 200 shown in Fig.1 in that an ink amount calculating means 111 
is used in place of the XYZ-RGB conversion means 106. 

[0039] Similarly to the color matching device 200 shown in Fig.1 , the color matching device 208 calculates the cor- 
30 responding color's X'. Y\ Z 9 in the reference side color image print device corresponding to the input color XYZ 5 from 
the original side color image display device 201 . 

[0090] If the corresponding color's X\ Y\ Z' 9 in the reference side color image print device 207 are calculated, the 
ink quantity calculating means 106 converts the corresponding color's X\ Y\ Z' 9 to the multi-color ink quantities to out- 
put the resulting ink quantities to the reference side color image print device 207. For converting tristimulus values XYZ 

35 to the multi-color ink quantities, the method described in (a) Japanese patent Application H-4-1 72246 (JP-A-6-22124) 
or (b) Japanese Patent Application H-6-79088 (JP-A-7-288706) may be used. However, other conversion methods may 
also be used. The entire disclosures of above (a), (b) are incorporated herein by reference thereto. 
[0091] With the color matching device or apparatus 204, 205 or 207, shown in Rgs.3, 4 and 6, in addition to the 
color matching device 200, it is possible to change the reproducing side color image display device 202 connected 

40 thereto to the reference side color image print device 207 by similar methods. 

[0092] Fig.8 shows a block diagram showing an illustrative structure of the hardware configuration for realizing the 
embodiments shown in Figs.1, 3, 5, 6 and 7. Referring to Fig.8, this hardware configuration is made up of a data 
processing device 400. including a storage (recording) medium 400, having the color matching program (program prod- 
uct) stored therein, an input device 401. an original side color image device 300 and a reproducing side color image 

45 device 301 . The storage medium 402 may be a magnetic disc, a semiconductor memory or the like memories. The 
present invention is carried out by the data processing device 400 reading out the color matching program for realizing 
the method shown in Fig.9 from the storage medium 402 and executing the read-out program. It is seHexplanatory that 
the storage medium 402 may be replaced by a communication network and the program product may be carried and/or 
supplied by a carrier medium such as a carrier wave having physical nature. 

so [0093] The meritorious effects of the present invention are summarized as follows. 

[0094] According to the present invention, as described above, since color matching is effected in meeting with the 
incomplete chromatic adaptation state in the human visual system under the conditions of different reference whites, it 
is possible to match the color appearance of a color image between color image devices having different reference 
whites, such as scanners, monitors or printers. 

55 [0095] Moreover, with the color matching method of the present invention, the color conversion can be realized at 
a high speed despite low computational cost and load. 

It should be noted that other objects and aspects of the present invention will become apparent in the entire disclosure 
and that modifications may be done without departing the gist and scope of the present invention as disclosed herein 
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and claimed as appended herewith. 

Also it should be noted that any combination of the disclosed and/or claimed elements, matters and/or items may fall 
under the modifications aforementioned. 

5 Claims 

1. A color matching method for converting an optional color in an originating color space to a color in a target color 
space, keeping the color appearance, characterized by comprising the steps of: 

10 (a) restoring spectral power distribution characteristics of an originating color space from a correlated color 

temperature of an originating reference white which is the reference white in the originating color space, and 
restoring spectral power distribution characteristics of a target color space from a correlated color temperature 
of a target reference white which is the reference white in said target color space; 

(b) restoring, as a surface reflectance of the originating reference white, a first white surface reflectance from 
75 tristimulus values of the originating reference white, said spectral power distribution characteristics of the orig- 
inating color space and human color matching functions, and also restoring a second white surface reflectance 
from tristimulus values of said originating reference white, spectral power distribution characteristics of said tar- 
get color space and the human color matching functions; 

(c) obtaining an adaptive white surface reflectance through interpolation of the first and second white surface 
20 reflectances; 

(d) calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance to obtain a 
spectral chromatic adaptation ratio; 

(e) restoring a surface reflectance of said optional color in said originating color space by employing tristimulus 
values of said optional color, spectral power distribution characteristics of said originating color space and the 

25 human color matching functions; 

(f) multiplying said surface reflectance of the optional color with said spectral chromatic adaptation ratio to 
obtain an adaptive surface reflectance of said optional color; and 

(g) finding tristimulus values of the color in said target color space from the adaptive surface reflectance of said 
optional color, the spectral power distribution characteristics of said target color space and the human color 

30 matching functions. 

2. The color matching method as defined in claim 1 

wherein in said step (c) of obtaining adaptation white surface reflectance through interpolation of said first 
and second white surface reflectances, an incomplete chromatic adaptation coefficient is changed responsive to 
35 the state of incomplete chromatic adaptation in the human visual system, the degree of interpolation of said first 
and second white surface reflectances is freely changed to obtain said adaptation white surface reflectance. 

3. The color matching method as defined in claim 2, further comprising the step of: calculating an optimum incomplete 
chromatic adaptation coefficient from combination of the reference white of said originating color space and the ref- 

40 erence white of said target color space, said incomplete chromatic adaptation coefficient obtained from said opti- 
mum incomplete chromatic adaptation coefficient calculating step being used. 

4. A color matching apparatus for converting an optional color in an originating color image device to a color in a target 
color image device, keeping the color appearance, characterized by comprising: 

45 

(a) means for restoring spectral power distribution characteristics of an originating color space from a corre- 
lated color temperature of an originating reference white which is the reference white of said originating color 
image device, and restoring spectral power distribution characteristics of the target color space from a corre- 
lated color temperature of a target reference white which is the reference white of said target color image 

so device; 

(b) means for restoring, as a surface reflectance of said originating reference white, a first white surface reflect- 
ance from tristimulus values of the originating reference white, said spectral power distribution characteristics 
of said originating color space and human color matching functions, and also restoring a second white surface 
reflectance from tristimulus values of said originating reference white, spectral power distribution characteris- 

55 tics of said target color space and human color matching functions; 

(c) means for obtaining an adaptive white surface reflectance through interpolation of the first and second 
white surface reflectances; 

(d) means for calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance 
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to obtain a spectral chromatic adaptation ratio; 

(e) means for restoring a surface reflectance of said optional color in sad originating color image device by 
employing tristimuJus values of said optional color, spectral power distribution characteristics of said originating 
color space and the human color matching functions; 
5 (f) means for multiplying said surface reflectance of the optional color with said spectral chromatic adaptation 

ratio to obtain an adaptive surface reflectance of said optional color; 

(g) means for finding tristimulus values of the color in said target color image device from the adaptive surface 
reflectance of said optional color, the spectral power distribution characteristics of said target color space and 
the human color matching functions, and 
w (h) means for converting the tristimulus values of the color in said target color image device into device color 

data of said target color image device. 

5. The color matching apparatus as defined in claim 4 

wherein in said means for obtaining adaptive white surface reflectance through interpolation of the first and 
is second white surface reflectances, an incomplete chromatic adaptation coefficient is changed responsive to the 
state of incomplete chromatic adaptation in the human visual system, the degree of interpolation of said first and 
second white surface reflectances is freely changed to obtain said adaptive white surface reflectance. 

The color matching apparatus as defined in claim 5, further comprising: 

means for calculating an optimum incomplete chromatic adaptation coefficient from combination of the refer- 
ence white of said originating color image device and the reference white of said target color image device, 
fcaid incomplete chromatic adaptation coefficient obtained from said optimum incomplete chromatic adaptation 
coefficient calculating means being used. 

The color matching apparatus as defined in any one of claims 4 to 6. 

wherein a color image display device is provided in the originating image device and a color image display 
device is provided in the target color image device. 

30 8. The color matching apparatus as defined in any one of claims 4 to 6. 

wherein a color image display device is provided in the originating image device and a color image print 
device is provided in the target color image device for color reproduction with multiple color inks; 

there being provided means for calculating ink quantity calculating as means for converting the tristimulus val- 
35 ues of the color in said target color image device. 

9. The color matching apparatus as defined in any one of claims 4 to 8, 

wherein there is provided a color image input device in the originating color image device by providing 
means for calculating the correlated color temperature of the reference white from chromaticity of the reference 
40 white of the color image device. 

10. A machine-readable program product performing, on a computer, the following program steps comprising: 

(a) restoring spectral power distribution characteristics of an originating color space from a correlated color 
45 temperature of an originating reference white which is the reference white of an originating color image device, 

and restoring spectral power distribution characteristics of a target color space from a correlated color temper- 
ature of a target reference white which is the reference white of a target color image device; (b) restoring, as a 
surface reflectance of said reference white, a first white surface reflectance from tristimulus values of the ref- 
erence white, said spectral power distribution characteristics of said originating color space and human color 
so matching functions, and also restoring a second white surface reflectance from tristimulus values of said orig- 

inating reference white, spectral power distribution characteristics of said target color space and human color 
matching functions; 

(c) obtaining an adaptive white surface reflectance through interpolation of the first and second white surface 
reflectances; 

55 (d) calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance to obtain a 

spectral chromatic adaptation ratio; 

(e) restoring a surface reflectance of said optional color in said originating color image device by enrploying 
said tristimulus values of said optional color, spectral power distribution characteristics of said originating color 
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space and the human color matching functions; 

(f) multiplying said surface reflectance of the optional color with said spectral chromatic adaptation ratio to 
obtain an adaptive surface reflectance of said optional color; 

(g) finding tristimulus values of the color in said target color image device from the adaptive surface reflectance 
5 of said optional color, the spectral power distribution characteristics of said target color space and the human 

color matching functions; and 

(h) converting the tristimulus values of the color in said target color image device into device color data of said 
target color image device. 

io 1 1. A machine-readable medium having said program product as defined in claim 10. 

12. A machine-readable program product performing, on a computer, 

(a) computing correlated color temperature of a reference white from chromaticity of a reference white of a 
15 color image device; 

(b) restoring spectral power distribution characteristics of an originating color space from a correlated color 
temperature of an originating reference white which is the reference white of an originating color image device, 
and restoring spectral power distribution characteristics of a target color space from a correlated color temper- 
ature of a target reference white which is the reference white of a target color image device; 

20 (c) restoring, as a surface reflectance of said reference white, a first white surface reflectance from tristimulus 

values of the reference white, said spectral power distribution characteristics of said originating color space 
and human color matching functions, and also restoring a second white surface reflectance from tristimulus 
values of said originating reference white, spectral power distribution characteristics of said target color space 
and human color matching functions; 

25 (d) obtaining an adaptive white surface reflectance through interpolation of the first and second white surface 

reflectances; 

(e) calculating a ratio of the first white surface reflectance to the adaptive white surface reflectance to obtain a 
spectral chromatic adaptation ratio; 

(f) restoring a surface reflectance of said optional color in said originating color image device by employing said 
30 tristimulus values of said optional color, spectral power distribution characteristics of said originating color 

space and the human color matching functions; 

(g) multiplying said surface reflectance of the optional color with said spectral chromatic adaptation ratio to 
obtain an adaptive surface reflectance of said optional color; 

(h) finding tristimulus values of the color in said target color image device from the adaptive surface reflectance 
35 of said optional color, the spectral power distribution characteristics of said target color image device and the 

human color matching functions; and (i) converting tristimulus values of the color in said target color image 
device into device color data of said target color image device. 

13. A machine-readable medium having said program product as defined in claim 12. 

40 
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